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Abstract 

We study the specific heat dependence of the periodic Anderson model (PAM) in the limit of 
U = oo employing the X-boson treatment in two different regimes of the PAM: the heavy fermion 
Kondo (HF-K) and the local magnetic moment regime (HF-LMM). We obtain a multiple peak 
structure for the specific heat in agreement with the experimental results as well as the increase 
of the electronic effective mass at low temperatures associated with the HF-K regime. The entropy 
per site at low T tends to zero in the HF-K regime, corresponding to a singlet ground state, and it 
tends to ksln2 in the HF-LMM, corresponding to a doublet ground state at each site. The linear 
coefficient 7(T) = C V (T)/T of the specific heat qualitatively agrees with the experimental results 
obtained for different materials in the two regimes considered here. 

PACS numbers: 71.10.-w, 71.27.+a, 71.28.+d, 75.20.Hr, 75.30.Mb, 75.40. Cx 
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I. INTRODUCTION 



Heavy fermions metals (HF) are intermetallic compounds of certain lanthanides and 
actinides. They contain a lattice of / ions, embedded in a Fermi sea of conduction electrons c. 
These materials exhibit very large specific heats and spin susceptibilities at low temperatures, 
and these characteristics are associated with high effective masses of quasi-particles, caused 
by the local strong Coulomb interaction between the / electrons. The periodic Anderson 
model (PAM) is considered a fairly good description of these systems, and the relevant 
parameters of this model are: the hybridization V between the localized /-electrons and the 
itinerant c electrons (associated with the s, p and d orbitals), the value of the energy Ef of 
the localized / level, the value U that determines the Coulomb repulsion energy between 
the / electrons, the dispersion relation of the c-electrons and the chemical potential \i. The 
physics of these materials is dominated by the competition between the Kondo effect, that 
tends to quench the magnetic moments of the / levels,— and the RKKY (Rudermann-Kittel- 
Kasuya-Yosida) interaction, that tends to introduce antiferromagnetic order in the latticed 
Some of these compounds develop superconducting states at low temperatures, and their 
nature is not yet completely understood.^. 

The thermodynamic properties of the PAM have been recently studied employing different 
methods: equation of motion (EOM) 6 ' 7 and two variants of the mean field theory^ that 
respectively describe the weak coupling limit (Hartre-Fock) and the strong coupling one 
(slave boson mean field theory (SBMFT) 9 ' 10 ! 11 ). An interesting result in all these works&^ 
is the presence of a multiple peak structure in the temperature dependence of the specific 
heat. In the works that use the EOM, the first peak was associated with spin fluctuations 
and the second one with charge fluctuations^. Bernhard et.al.—, defined a lattice Kondo 
temperature Tk (we call it Tkl i n this work) as equal to the T corresponding to the minimum 
of the temperature derivative of the average (cl a fi a \, that measures the "transference" of 
electrons from the localized levels to the conduction band and vice- versa. On the other hand, 
Peres et.al.—, define a correlation temperature T* (we call it To in this work) equal to the 
T corresponding to the maximum of the temperature derivative of the parameter (c^fia)- 



As they employ the SBMFT the dynamical fluctuations do not appear explicitly in their 
Green's Functions (GF), and they concluded that the multiple peak structure of the specific 
heat is only associated with band shape effects. We believe that there is no antagonism 
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between these two explanations: the dynamical fluctuations modify the bands (even in the 
SBMFT through their mean field contributions to the GF), and it is their spectral densities 
near the Fermi surface that determine the specific heat. 

The dynamical mean field theory*^. (DMFT), with the help of the numerical renormaliza- 
tion group (NRG), has been recently employed to study the PAM^ 3 - and the Kondo lattice.— 
These works establish the simultaneous existence of two scales of energy: one high temper- 
ature scale associated with conventional Kondo scattering and a lower temperature scale, 
usually called coherence temperature, that is associated with the formation of quasi-particles 
with strongly enhanced effective masses.— The same result has been also obtained by ap- 
plying the SBMFT in the large-TV approach to the Kondo lattice.— The notation of these 
temperature scales is not uniform in the literature, and we shall then use T for the lattice 
coherence temperature , 13 ! 14 and introduce Tkl for the high temperature scale, leaving Tk 
for the impurity Kondo temperature. 

In the present work we study the thermodynamic properties of the PAM in the limit 
of U — ► oo employing the X-boson metho d 1 ^ 1 ^ 1 ^ two regimes of the PAM: the heavy 
fermion Kondo regime (HF-K) and the local magnetic moment regime (HF-LMM).- 4 ^ The 
spurious second order phase transitions of the slave-boson approach , 16 > 19 at intermediate 
temperatures and when Ef « //, are not present in the X-boson method. This method 
also includes explicitly some temporal quantum fluctuations in the GF,~ and we analyze 
some possible consequences of this difference. We show that the parameters T and T KL 
are determined by the position of the maximum of Cy vs. T, and we discuss the physical 
origin of these maxima. We show that at low T the entropy per site tends to zero in the 
HF-K regime, corresponding to a singlet ground state, while it tends to ksln2 in the HF- 
LMM, corresponding to a doublet ground state at each site and opening the way for the 
possible formation of an antiferromagnetic order (AF) at low temperatures (T < To). We 
show that when the system is in the HF-K regime and the system goes at low temperatures 
through the Kondo resonance, there is a huge increase of the effective mass of the f electrons, 
associated with the linear coefficient 7 Q of the specific heat. Our calculations of 7(T) = C v /T 
qualitatively agree with the experimental results at low temperatures for some materials 
classified in the two regimes considered. 
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II. MODEL AND METHOD 



To study the PAM Hamiltonian in the U = oo limit, the f-levels with double occupa- 
tion are projected out from the space of the local states by employing the Hubbard X 
operators , 20 i 21 i 22 and we obtain 



k,cr a,f 

+ E (^^4o/V + v f \A,« x f,°°) ■ (!) 

k,cr 

The first term of the equation represents the Hamiltonian of the conduction electrons (c- 
electrons), associated with the itinerant electrons (s, p and d orbitals). The second term 
describes the localized / levels and the last one corresponds to the interaction between the c- 
electrons and the /-electrons via hybridization between the / and c states. This Hamiltonian 
can be treated by the X-boson approac h 16 ' 17 ' 18 for the lattice case, and the cumulant Green's 
function (GF) are then given by 

Gg{z) = , „ - D ° , (2) 

-E fj<7 ) ( z - £]l(7 ) - \V a (k)\W a 

= 7 , (3) 



= 7 = r DMk) , (4) 



and 

G&(z) = - 

z-E f ^ ( z - £]l(7 ) - \V a (k)\W a 
where z = to + ir], with 77 — ^ + . The correlations appear in the chain X-boson approach 
through the quantity D a = R + nf >a , where R = (X 0t0 ). This simple factor introduces 
essential differences with the MFSB T 16 ' 17 and these GF can not be transformed into those 
of two hybridized bands of uncorrelated electrons with renormalized parameters, as it is 
done in the MFSBT. A renormalized Ef = Ef + A also appears here in the X-boson method, 
with A given by 

S k,<7 

n F Ok,<r(+)) - n F (a; kj(T (-)) 



A 

k.i 



(e Ka -E f ^ + 4|K(k)| 2 L> (J 



(5) 
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where rip (x) is the Fermi-Dirac distribution and N s is the number of sites. To simplify the 
calculations we shall consider a rectangular conduction band of width W = 2D, centered at 
the origin, and a real hybridization constant V^(k) = V. We then obtain 

where the values c<j k CT (±) are the poles of the GF, given by 



^ k>CT (±) = - (e k ^ + Ef 



±y(e k , a -E f y + 4V*D a . (7) 

In the X-boson approach D a = R + rif a must be calculated self-consistently through the 
minimization of the corresponding thermodynamic potential with respect to the parameter 
R . When the total number of electrons N to t, the temperature T and the volume V s are 
kept constant one should minimize the Helmholtz free energy F, but the same minimum 
is obtained by employing the thermodynamic potential Q = —kBTln(Q), (where Q is the 
Grand Partition Function) while keeping T, V,, and the chemical potential fi constant (this 
result is easily obtained by employing standard thermodynamic techniques). In the X-boson 
method the grand thermodynamic potential i o iei17 

n = n + ^ ln t 1 + e M-P ^k, CT (±))] 

k,o-,± 

+ N a A(R-l), (8) 

where 



fi = --ln 



l + 2exp(-pE f ) 



(9) 



[1 + eM-PEf)) 2 

(3 = 1/ksT, and ks is the Boltzmann constant. By standard thermodynamic techniques we 
find 

S=-(%) , (10 ) 



and 



\dT J N v \dTj. 



c " ~ 1 ' or) ~ 1 \lm , 

U1 / N tot ,Vs V U1 / ^V s 



+ T\g) f^l . (11) 



dT J N v V 8T 



fl,Vs 
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For the rectangular conduction band we find from Eqs. ()8l9j) that, in the absence of magnetic 
field, 



where 



and 



+ 



T 
D 



d 2 Q 

~&T 2 , , 

2 2 I'D 



-T 



d 2 n 
"df 2 



y / dxu 2 {x) np (u>g(x)) [1 — np {ui{x))\ 



T N„ 



d\A(R-l)) 
dT 2 



w±{x) = -[x + Ef)±-\[x-Ef) +A\VrD a 



d 2 n 

dT 2 



-2N S k B (3 2 Ej exp(/3E f ) 



[3 + 2exp(/?E / )] 



2 ' 



exp(pE f ) + 2 exp(pE f ) + 1 
We have calculated T (d[i/dT) NtgtVs (dN tat /dT) Vs numerically. 



(12) 



(13) 



(14) 



III. RESULTS AND DISCUSSION 

A schematic classification for typical Ce-based compounds was given earlier by Varma 23 
and reintroduced by Steglich et. al4«^ in terms of the dimensionless coupling constant for 
the exchange between the local / spin and the conduction-electron spins, g = Np \J\, where 
Np is the conduction-band density of states at the Fermi energy and J is connected to the 
parameters of the PAM via a Schrieffer- Wolff transformation,— J = \2V 2 / ( E f — (for 
U — > oo). The behavior of these compounds can then be qualitatively organized through 
the parameter 

g = 2V 2 p c (fi)/\E f -fi\, (15) 

where p c (/i) = 1/2D is independent of \i in our case, Ef is the bare localized / level and 
we employ the chemical potential fi obtained in the self-consistent X-boson calculation for 
the total number of particles per site n tot = 2{rif a + n ca ) at the very low temperature 
T = 0.000LD. 
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When g > 1, the compound under consideration is in the intermediate valence (IV) region, 
while for g < 1 it is in the heavy fermion Kondo regime (HF-K). There exists a critical value 
g c at which the Kondo and the RKKY interactions have the same strength, and non Fermi- 
liquid (NFL) effects have been postulated for systems with g = g c . For g c < g < 1 the 
magnetic local moments are quenched at very low temperatures and the system presents a 
Fermi liquid behavior, while for g < g c we have the local magnetic moment regime (HF- 
LMM). We point out that the parameter g classifies the regimes of the PAM only in a 
qualitative way. Finally, as in its present form the X-boson approach includes hybridization 
effects only to second order in V and the self-energy does not depend on the wave vector,— 
the RKKY effects are not present and we cannot discuss non Fermi-liquid behavior within 
the present approximation, nor find the value of g c . 

A. Results for the HF- regime 
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FIG. 1: Specific heat Cy and occupation of the localized / levels nj vs. temperature T, for the 
HF-K regime. Notice that the intensity of the two maxima of Cy at lower T first decreases and 
then disappears when the number n c = 1.6 — rif of conduction electrons increases, while their 
positions To and Tkl shift to the right. 

In Fig. [I\we show the specific heat C v and the occupation rif of the localized level as a 
function of the temperature T for several sets of parameters that describe a HF-Kondo situa- 
tion (HF-K): hybridization V = 0.2D, total number of particles per lattice site n tot = 1.6 and 
bandwidth W = 2D; the bare localized energy levels Ef employed and the corresponding 
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n m =1 .6, V=0.2D, W=2D 




Steglich's parameters'^ g (for T = 0.0001D) are E f = -0.05D, g = 0.426; E f = -0.10D, 
g = 0.392; E f = -0.15D, g = 0.365 and E s = -0.20D, g = 0.342. In Fig. Hwe plot both 
the / and c spectral densities for Ef = — 0.15.D, showing the high density of states / at the 
chemical potential for that set of parameters. 



FIG. 2: The / and c electronic density of states pf jCr (u}) and p c ^{oj) vs. u> at low temperatures, for 
the HF-K regime. The position of the renormalized energy level E = Ef + A and of the chemical 
potential p are shown. Similar results were obtained for other set of parameters in the same regime. 

In Fig. |3] we plot the specific heat, the entropy and the total occupation of the / level as 
a function of T for the same set of parameters of Fig. |21 We can see that the entropy S per 
site tends to zero at low temperatures, and we conclude that the system goes to the Kondo 
singlet ground-state, while the huge value of density of states / at the chemical potential \i 
(see figure |2J) signals the simultaneous appearance of the Kondo resonance. The inset shows, 
in agreement with reference that an entropy approximately equal to n c ks ln(2) is lost 
below T KL . 

At high temperatures the system behaves like a collection of independent / and c elec- 
trons, because both V/ksT and D/ksT tend to zero like 1/T. For a system with a chemical 
potential that remains finite when T — ► oo the entropy per site would tend to ks ln(12), 
corresponding to the dimensionality of the purely local Hilbert space, but this result is not 
true any more when we keep fixed the total number of electrons per site n to t, as discussed in 
more detail in Appendix El The value of the entropy when ntot = 1.6 (the value employed 
in figures Q]|ni) is equal to 0.998765 ks lnY2 and is indistinguishable from ks /n(12) in the 
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figures. 



FIG. 3: Entropy S, Specific heat C v and total occupation of the localized levels nj vs. T, for 
the same set of parameters of Fig. |2I The behavior of the curves is similar for other parameters 
corresponding to the same HF-K regime. 

In the Fig. 0] we plot the specific heat C v , the chemical potential /i, the temperature 
derivative of the parameter (c\ a fi a y = rif Cta , the occupation number rif of the / levels, and 
the value at /i of the / density of states pf^(/i), as a function of T for a system in the HF-K 
regime. The average i^iafia^j = nj c ^ a is associated with the GF G^ a ( y z n ) ) and it is a measure 
of the mixing of the c and / states in the system. We associate the temperatures T and 
Tkl respectively to the maximum and the minimum of the (dnf c ^ a /dT) (they are named T* 
in Peres et.alA and Tk in Bernhard et.al^ respectively). They are both shown in the figure, 
where it is clear that T coincides with the first peak of Cy and Tkl with the second one. 
Only one extreme was observed in previous works.^ 

The chemical potential \x changes with temperature because we keep the total number of 
particles constant, and in our calculations it shifts into the hybridization gap at intermediate 
temperatures. The density of states at the chemical potential has a very large value 

at low T, but at T ~ To this quantity shows a drastic change with increasing T, taking a 
value close to zero when \i enters the gap, as shown in Fig. 0] The moderate increase in the 
average electronic energy that occurs in this process is the origin of the first peak of C v . 

It seems reasonable to use T as a measure of the Kondo lattice temperature, because 
crossing Tq (with decreasing T) we find that p/ )(T (/x) changes from very small to very large 
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values, and at the same time the entropy takes small values, tending to zero when T — > 
(cf. Fig. EI). These two properties characterize respectively the appearance of the Kondo 
peak and the formation of a singlet by the quenching of the magnetic moment of the / 
electrons by the cloud of c electrons, and our choice of T as the crossover parameter seems 
reasonable. 




FIG. 4: Specific heat C v , chemical potential //, density of states p/ )Cr (/x) of the /-electrons at ji, 
(dnf c ^ a /dT), and nj vs. T, for the same set of parameters given in Fig. ^but with Ef = — 0.15D. 
We indicate the temperatures Tj£ and Tkl*^ The inset shows the little peak in C v at T = To. The 
behavior for other sets of parameters corresponding to the same regime is similar. 

The temperatures To and Tkl coincide in our treatment with the first two peaks of 
CV, and correspond to the two energy scales that appear in the DMFT treatment of both 
the PAM^ and of the Kondo lattice.— The lattice coherence scale is associated in the 
first reference^ to the inverse effective mass at the Fermi surface, which in our treatment 
is proportional to the corresponding density of states: as shown in Fig. 0] this quantity 
develops a very large increase when T decreases below T . This same scale is also associated 
in the second reference^ to the onset of Fermi-liquid coherence, and can be determined 
from various dynamical quantities including the /-electron Kondo resonance.— It should 
be remarked that the drastic change that p/ CT (/i) suffers when T crosses T (cf. Fig. EJ), 
is consistent with the result reported in Section 3.5 of reference,— that the position of the 
peak of the resonance, "relative to the chemical potential shifts to below the Fermi level on 
a scale < To with increasing temperature." 

The second peak at T = T KL corresponds to a second and larger maximum of 
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corresponding to the maximum rate of change with T of the mixing be- 
tween / and c electrons. Even if the GF do not have explicit dynamic fluctuations, the 



changes in 



c ] f- 



measure a mean-field average of the charge fluctuations, and the peak 



is due to the energy necessary to change the average occupation of the different types of 

d(c\j la )/dT 



electrons. The crossover parameter T KL has to be used when d {c^f^j /dT shows a single 
extreme, as it happens in the half filled symmetric PAM^ 6 because /i is then fixed in the 
middle of the gap and the changes in population that originated T do not occur. 

The temperature Tjcl corresponds to the second energy scale in both DMFT works , 13 i 14 
and its corresponding energy is related to the transitions between the states of the pseudo 
gap that usually appears in that treatment. Our method only gives a real self energy, and we 
always obtain a real gap rather than a pseudo gap, but the meaning of Tkl remains the same. 
The energy associated to T KL has essentially the same origin as the one corresponding to 
the Kondo impurity Tk, describing the conventional incoherent Kondo scattering at higher 
temperatures.— In agreement with the results of several authors , 13 i 14 i 15 Tkl decreases when 
the number n c of conduction electrons decreases. 

Two scales are also found in the large- iV SBMFT of the Kondo lattice,^ characterized 
by the two maxima of Cy vs. T (though the second maximum is replaced in their Fig. 3b 
by a discontinuity at the phase transition that appears in the SBMFT at Tkl)- Although 
both the model and the treatment are different, the Cy maxima at T in our Fig. ^have 



the same trend they show in reference 15, namely that they decrease and finally disappear 
when the number n c of conduction electrons increases. 

Finally the third peak of Cy is located in a region corresponding to large T, and it 
measures the energy necessary to take the electrons from the low T distribution into the high 
temperature limit distribution, that depends on n tot (it would be the uniform distribution 
for n tot = 5/3, cf. Appendix |XJ) . As could be expected, the position of the high temperature 
peak is independent of the value of Ef, as it is shown in figure H 

In Fig. El we plot the coefficient j(T) = Cy(T)/T vs. T for several values of Ef corre- 
sponding to the HF-K regime. At low T the coefficient j(T) is proportional to the density 
of states at the chemical potential /i, and in simple cases this is proportional to the ef- 
fective mass m* of the associated quasi-particles. We see in that figure that the value 
of j(T)\t=o.o — lo at very low temperatures has a large increase when Ef changes from 
Ef = — 0.05.D to Ef = — 0.20-D, while the Steglich's parameter respectively decreases from 
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FIG. 5: Sommerfeld's coefficient j(T) = C V {T)/T vs. T at very low temperatures for the HF-K 
regime. We have 7(T')|t=o.o = To oc m*, where m* is the effective mass of the quasi-particles. The 
two vertical lines delimit the interval of values of To corresponding to the four different curves. In 
the inset one sees an anomaly in j(T) at T = Tq and verifies that T 7 << To for all the values of 
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FIG. 6: Sommerfeld's coefficient j(T) = C V (T)/T vs. T for the HF-K regime, at very low temper- 
atures (T ~ T 7 ). 

g = 0.426 to g — 0.342. This shows that the quasi-particles become heavier as the system 
evolves in the direction of the Kondo regime. The maximum of 7(T) is reached at a value 
T 7 that satisfies T 7 << T . The inset of Fig. shows j(T) vs. T/T for T/T ~ 1, and 
indicates an evident anomaly in 7(T) at T = Tq for all the values of Ef in the HF-K param- 
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eter region. The general form of our results agrees with the experimental measurements for 
HF systems,— i 27 i 28 i 29 i 30 i 31 i 32 as well as with those of a recent theoretical work.— 

In Fig. |U1 we present the same results of Fig. |5] for very low temperatures showing 
the maximum of 7(T) in detail. The results are qualitatively similar to the experimental 
measurements for different materials,— like C eAl ^ 28 ^ 9 ^ 30 ^ 2 and CePtSiA^ It is important to 
emphasize that the best qualitative agreement with the experimental results is given by the 
set of parameters that corresponds to the highest values of pf j0 -(fi) at very low temperatures, 
pointing out the strong HF character of these compounds. 

B. Results for the HF-LMM regime 

Now we consider several sets of parameters associated with the heavy Fermion local 
magnetic regime (HF-LMM), defined as before through the Steglich's criterion g < g c < 1^ 
but characterized in our treatment by its properties.- 2 ^ Our calculations in this regime show 
that the density of / states at the chemical potential p/ j(J (/x) for T —>■ has low values, that 
the specific heat has a two peak structure, and that the occupation number nf of localized 
states is close to 1. In Fig. [7|we plot the same properties given in Fig. ^for the following 
set of parameters, that correspond to the HF-LMM regime: hybridization V = 0.2D, total 
number of particles per site n tot = 1.9, bandwidth W = 2D, and localized energy level given 
by Ej = —0.20/}, Ef = —0.25D, and Ef = — 0.30D, corresponding to Steglich's parameters 
g = 0.216, g = 0.191, and g = 0.164 respectively. 

In Fig. [HI we plot both the / and c spectral densities for Ef = —0.20D and T = 0.001-D, 
showing the rather small value of the spectral density p/ j(7 (/i) at the chemical potential \i. 
Notice that in this case (HF-LMM regime) the renormalized energy Ef = Ef + A is well 
below the chemical potential /i. 

The Fig. El is similar to Fig. El but for the system in the HF-LMM regime. At very low 
temperatures the entropy S per site tends to ksln(2), pointing to the transformation of 
the singlet of the HF-K regime into a ground state consisting of a doublet at each site, that 
could be attributed to a spin 1/2 at each site. Like in the HF-K regime, the entropy S goes 
to a high temperature limit that depends only on n to t, as discussed in Appendix lA"l 

In Fig. El we present the same results of Fig. E]for the HF-LMM regime. In this regime 
the pf <a (fi) goes to zero only in the limit of very high temperatures, because /i is always in 
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FIG. 7: Specific heat C„ and occupation on the localized / levels nj vs. temperature T, for the 
HF-LMM regime. 
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FIG. 8: Density of states pf >a (u) and p Cia {uo) of / and c electrons as a function of u at low tem- 
perature, for parameters corresponding to the HF-LMM regime. The position of the renormalized 
energy level E = Et + A and of the chemical potential fi are shown. Other sets of parameters in 
the same regime show a similar behavior. 

the tail of the density of states (see figure IHJ) and does not cross through the gap, as it does 
in Fig. |U As a consequence, the peak corresponding to the one at the lowest T in figures [U 
El and |U is absent in this regime, and the first specific heat peak that appears is physically 
equivalent to the one at T^l m the HF-K regime, i.e.: it corresponds to a redistribution 
of states into higher energies. This Cy peak is less pronounced than in the HF-K region, 
because many of the low energy states are already occupied at small T because of the larger 
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value of N t . The second peak of Cy is now associated with the high temperature limit, and 
the same arguments employed in the HF-K region apply here; in particular its position is 
independent of the value of Ef (cf. Fig. [7J). 
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FIG. 9: Entropy S, specific heat C v and charge at the localized levels nj vs. T for a set of 
parameters corresponding to the HF-LMM regime. Other sets of parameters in the same regime 
show a similar behavior. 
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FIG. 10: Specific heat Cy, density of local states pf ;fT ((i) on the chemical potential fj,, drif c /dT and 
n* vs. T, for a set of parameters corresponding to the HF-LMM regime. Only the maximum of 
dnf c /dT coincides with a maximum of Cy at Tkl- Other sets of parameters in the same regime 
show a similar behavior. 

In Fig. El we present the same results of Fig. [7| (HF-LMM regime) but with the temper- 
ature scaled by Tkl- All the first peaks collapse at T = Tkl by construction, and the inset 
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shows a universal cubic dependence of C v vs. T at low temperatures. The same dependence 
was obtained by Tesanovic and Vails 3 - 3 by starting from the slave boson mean field treat- 
ment and considering the interaction between heavy fermions, mediated by fluctuations of 
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the X-boson method 



the slave bosons. As discussed after equation (All) in reference 
retains some time dynamics that is absent in the slave boson method, and this might lead 
to the same effect obtained by Tesanovic and Vails. 

Kim and Stewart^ conclude from their experimental results that, in some heavy Fermion 
compounds like CeCug, a significant contribution to the Cy comes from magnetic correla- 
tions. We are not explicitly including magnetic interactions into the model Hamiltonian, 
and as discussed before, our treatment neither includes the RKKY interaction. Therefore, 
we can not attribute this T 3 dependence of Cy to intersite magnetic interactions, but it 
would correspond to Sommerfeld's second tern in the Cy temperature expansion. 
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FIG. 11: Specific heat and occupation numbers of / electrons vs. T/Tkl in the HF-LMM regime. 
The inset presents the T/Tkl < 1 region, showing a cubic contribution C v oc j (T/Tjcl) + 
[i(T/T KL f. 

In Fig. ^]we show the coefficient 7(T) = C v /T vs. T for the HF-LMM regime, and in 
the inset we show 7(T) = C V {T)/T vs T/Tkl- The reduction of the effective masses m* oc 7 D 
corresponds to the crossover from the HF-LMM regime to Steglich's stable moment regime 
(SMR ) corresponding 4 ' 5 to g << 1. The reduction of the effective mass is again associated 
to the decrease of Pf, a {n) at very low temperatures (T < T KL ). 

In Fig. EI] we present in a linear scale the same results of Fig. ^]at T ~ Tkl, but only 
for Et = —0.20D. The inset shows the linear behavior of (j(T) — 7 G ) vs. T 2 , corresponding 
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FIG. 12: Sommerfeld's coefficient 7(T) = C V {T)/T vs. T in the HF-LMM regime at very low 
temperatures. The coefficient 7(T)|t=o.o = 7o oc m*, where m* is the reduced mass of the quasi- 
particles. 



FIG. 13: Sommerfeld's coefficient 7(T) = C V (T)/T vs. T in the HF-LMM regime at very low 
temperatures. The inset shows the linear behavior of the (7(T) — 7 D ) vs. T 2 associated with 
Sommerfeld's second term in the T expansion of Cy . 

to the cubic dependence of Cy discussed above. 
1. Discussion 

It is clear that our method does not give the Kondo peak nor the quenching of the 
magnetic moment in the HF-LMM region. This is the region where the RKKY interactions 
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override the Kondo effect, and one is not quite clear of what should be the exact paramagnetic 
solution in that region. The Kondo resonance is observed in some theoretical calculations 
of the paramagnetic region with a moderately small value of Steglich's parameter g 4 ' 5 , like 
in Fig. 3f of reference {?;, but we are not certain that these calculations correspond to 
the HF-LMM region. We must also consider that we have not included explicit magnetic 
interactions into the Hamiltonian, and also that there is no RKKY interaction because 
only second order cumulants are employed in the basic CHA expansion, as shown by the 
corresponding diagrams (see e.g. Fig. 1 or reference 1 16T). 

The absence of fourth and higher order cumulants in the CHA leads to a disturbing 
question: how do we obtain any Kondo resonance at all in the HF-K region, if apparently 
there are no spin flip processes in our calculation? To answer this question we notice that 
although the X-boson method starts from an expansion without spin flips, it is followed by 
a minimization of the free energy with respect to R = (X 00 ), with the completeness as a 
subsidiary condition (cf. Eq. (4) of reference [3) . In this procedure the original solution is 
modified so that it gives the best possible physical solution (as regards to the free energy 
F) compatible with the general form of the CHA. The shifting of the localized energy Ef 
to the renormalized value Ef = Ef + A optimizes the free energy F giving at the same time 
a peak near the chemical potential /i, that takes the place of the Kondo resonance when 
parameters in the HF-K region are used. This result is not achieved in the HF-LMM region 
because the required shift Ef would not minimize F, and one is then left with an unquenched 
doublet without the corresponding spectral density increase around /i. A complementary 
explanation is that the constrained minimization would correspond to adding to the CHA 
expansion an extra number of undefined higher order diagrams, that would result in the 
formation of the Kondo resonance and of the quenching of the magnetic moment by the 
conduction electrons in the HF-K region. 

We should point out that the same type of arguments presented here would also apply to 
the SBMFT, because after making the mean field approximation of the slave boson operator, 
one is left with two hybridized bands without any explicit correlation, and it is again the 
free energy restrained minimization that leads to the formation of a resonance that describes 
the Kondo peak. 
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IV. SUMMARY AND CONCLUSIONS 



We employed the X-boson technique to calculate the specific heat of the PAM in two 
regimes (HF-K and HF-LMM), and we find agreement with recent theoretical wor k 6 i 7 i 8 i 13 i 14 i 15 
as well as qualitative agreement with measurements^ of Sommerfeld's coefficient j(T) = 
C V (T)/T. The low temperature entropy shows a crossover from a singlet ground state in the 
HF-K regime to a lattice of doublets in the HF-LMM region, and we discuss the origin of 
this behavior. We conclude that our method gives fairly good results in the Kondo region 
not very far from the intermediate valence region, as well as in this last region. 

The specific heat in the HF-LMM region has a C V (T) ~ /3T 3 contribution, that can 
be attributed to Sommerfeld's second term in the temperature expansion of Cy. It would 
be interesting to calculate the electrical resistivity in the two regimes at low temperatures 
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APPENDIX A: THE ENTROPY OF THE PAM FOR VERY HIGH T 

When T —>■ oo one is tempted to consider a uniform occupation of all the states, but 
we shall show that this is only true when the chemical potential is kept constant. This 
assumption is not generally true for a system with a fixed total number n to t = N tot / N$ of 
electrons per site (iVg is the number of sites) , and we shall calculate the entropy in this limit 
as a function of n to t ■ 

When T — > oo we can set (we take ks = 1) V/T = W/T = 0, and we can write the 
thermodynamic potential per site 



(T<T ). 



T In [(1 + 2 exp [- (E f - /T]) 



x (l + exp[+/i/T]) 2 ] 



(Al) 
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FIG. 14: Entropy per site of the PAM for T — > oo as a function of the total number of electrons 
per site y = n to t- The horizontal line corresponds to St-»oo (5/3) = In (12). 



and 



For convenience we write 



2 exp [- (E f - m) /T] 2 exp [+/i/T] 
"" " 1 + 2 exp [- - M ) /T] 1 + exp [+///T] ' 



x = exp [— /x/T] , 
y = n tot , 

e f = exp [Ef/T] (A3) 



so that Eq. (jA2j) can be written as 



2 2 

S/= + — =-. (A4) 

y e f x + 2 x + 1 K ' 

As we are interested in the limit T — > oo we write here e/ = 1, an solve x as a function of y: 



(4 - 3y + yi6" 



y 2 



x = ^ ^. (A5) 



The solution with — a/16 + |/ 2 gives negative values of x and we have to bound the total 
number of electrons n tot — y to the interval [0, 3], i.e. only take 

< y < 3. (A6) 
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It is clear that except for x = 1 (which corresponds to n tot = 5/3) we have // = a T, 
that goes to infinity together with T, and from Eq. (jA4|) or Eq. (jA5|) follows that the value 
of a has the sign of {n tot — 5/3). 

The thermodynamic potential and the entropy can now be calculated for T — ► oo as a 
function of y = n tot . The expression for the entropy is 



s. 



T^oo 



(y) = y ln 



4-3y+ y/16 + y 2 
2y 

2 



(A7) 



(A-y+ ^16 + y 2 ) (4 + y + ^16 + y 

(4 - 3y + Vie + r 



+ In 



and is plotted in figure ITU 

The entropy for n tot = 1.6 is <5r->-oo (1-6) = 0.998765 ln (12) (see figureEJ) while for n tot = 1.9 
is SY^oc (1.9) = 0.984597 ln (12) (see figure^. It follows from Eq. (jSU) that for T -> 00 it 
is 



where as before x = exp [— fi/T]. 



2 + x 



and n r 



1 + 



(A8) 
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